Journal of Systematics and Evolution

Comparative phylogeography provides a powerful framework to identify the historical processes shaping biodiversity hotspots by testing whether co-distributed species exhibit shared patterns of diversification. Southeast Asia harbors exceptional biodiversity, yet the extent to which common paleogeographic and climatic drivers have structured diversification across taxa remains poorly understood. Here, we investigated the evolutionary history of five widespread Microhyla species complexes distributed across the Indochinese Peninsula and adjacent regions using dense mitochondrial sampling (1,388 individuals) combined with genome-scale ddRAD sequencing (280 individuals). Across all complexes, both approaches recovered multiple lineages and remarkably congruent phylogeographic breaks and contact/intergradation zones associated with major Indochinese regions, including Myanmar, the Tenasserim-Malay Peninsula, southern Vietnam, and northern Vietnam-southern China, supporting the hypothesis that Indochina functions as a mosaic of stable biogeographic units. However, lineage divergence varied substantially among complexes, suggesting that common biogeographic drivers interacted with species-specific demographic histories. Phylogenomic analyses and cyto-nuclear discordance further revealed historical introgression in every complex, indicating that diversification involved both long-term allopatric isolation and reticulate evolution. These findings portray the Indochinese biodiversity hotspot as a dynamic evolutionary system where cycles of fragmentation and reconnection have repeatedly reshaped lineage boundaries. Moreover, the complex phylogeographic structure recovered exemplifies the urgent need for taxonomic revisions, for which genome-scale data provide an essential framework to validate mitochondrial hypotheses and assess admixture patterns for species delimitation. Finally, regions such as the southern Annamites and Tenasserim Hills emerged as recurrent hotspots of genetic diversity across independent lineages, highlighting their importance for conserving not only species/lineage richness but also the evolutionary processes and adaptive potential that sustain biodiversity.

Biodiversity aggregators such as GBIF provide unprecedented access to global biodiversity data, yet their representativeness remains uneven across space and taxa. This study examined the spatial and taxonomic structure of global vascular plant data available on GBIF. Six filters were applied to the GBIF vascular plant dataset, resulting in the removal of 54% of all records. Together, the filters explained more than 90% of the identified spatial issues, with duplicate and missing coordinates accounting for most of the variation. A higher number of occurrence records was associated with a greater number of spatial issues. Record distributions became progressively more even at finer taxonomic levels, from orders to species. The time series of occurrences for species, genera, and families increased sharply after 1800 and continued to rise, with no apparent stabilisation. Of the 824 ecoregions covered, 73 accounted for 72% of all occurrence records. These ecoregions spanned all continents but were strongly concentrated in Europe, followed by North America and Oceania. The analyses reveal four key patterns: (1) data volume is positively associated with spatial issues; (2) a small number of taxa account for a large proportion of records, whereas many are represented by relatively few; (3) occurrence data aggregated by GBIF have increased continuously since 1800; and (4) record coverage remains highly uneven across the world's ecoregions. These results highlight the substantial contribution of biodiversity data aggregators to expanding access to biological information while demonstrating the persistent spatial and taxonomic biases that shape their contents. Such biases should be explicitly considered when assessing data completeness and quality and when using aggregated occurrence records to infer global biodiversity patterns.

Determining the correct donor species/lineages of grass-to-grass lateral gene transfer (LGT) is vital for deducing specific donor features that could help inform the mechanism of transfer. This requires a dataset spanning a broad range of species to achieve the phylogenetic resolution necessary for precise donor inference. As grass-to-grass LGT often involves the transfer of multi-gene DNA fragments, they can contain additional sequences that allow for accurate orthologous comparisons, such as nuclear DNA of plastid origin (NUPTs). Here we systematically scan for NUPTs in the genomes of four Alloteropsis semialata accessions, whose LGTs have previously been characterised. Using the abundant Panicoideae chloroplast sequences, we reconstruct NUPT phylogenies and infer two lateral acquisitions: one from Paniceae/Digitaria and another from Andropogoneae/Eremochloa adjacent to a previously identified LGT. We then assembled and included an additional 12 Eremochloa chloroplast genomes in the analysis and showed the likely donor was Eremochloa attenuata. Subsequent short-read mapping from E. attenuata to the nuclear region flanking this NUPT showed consistent coverage across the region, including the previously identified LGT, supporting co-transfer. Overall this study highlights the potential for NUPTs to better identify the donors of grass-to-grass LGT.

Sinia rhodoleuca, the sole species of the monotypic genus Sinia (Ochnaceae), was previously transferred to Sauvagesia based mainly on morphological similarities. However, its phylogenetic position has remained unresolved because molecular data for the species were unavailable. Here, we generated genomic data for Sinia rhodoleuca and reconstructed its phylogenetic position within Sauvagesieae. Our phylogenomic analyses consistently recovered Sinia rhodoleuca as sister to Indosinia, whereas the Neotropical Sauvagesia formed a distantly related lineage, rendering Sauvagesia broadly circumscribed non-monophyletic. Comparative morphological evidence further supports the close relationship between Sinia and Indosinia, particularly in their closely parallel secondary veins, lacerate stipules, and prominent petaloid staminodes, while differences in floral characters support their recognition as distinct genera. We therefore reinstate Sinia as a distinct genus and provide a revised taxonomic treatment of Sinia rhodoleuca. Our study demonstrates how phylogenomic evidence can resolve long-standing taxonomic uncertainties and reveal evolutionary relationships obscured by morphological similarity.

Telomeres are nucleoprotein structures that protect chromosome ends and are maintained by the Telomerase Reverse Transcriptase (TERT) protein that uses a noncoding Telomerase RNA (TR) as a template. In monkeyflowers, Mimulus lewisii had an ancient TR gene duplication, synthesizing an evolutionarily atypical sequence heterogeneous telomere. How TERT interacts with both TR paralogs during telomere maintenance is unknown and answers can shed novel insights underlying telomere function. Using new genome assemblies we discovered TERT is rapidly evolving in lineages sharing the TR duplication. We investigated the functional consequences arising from the rapid evolution, first by using yeast three-hybrid and testing the physical binding between conspecific and heterospecific TERT-TR combinations. Results showed TERT binds both ancestral (TR1) and derived (TR2) TR paralogs in M. lewisii, but not in species without a functioning TR2. We located the region of TR binding to amino acids near the KRxR motif. We then combined next-generation sequencing with Telomeric Repeat Amplification Protocol and discovered M. lewisii had high telomerase activity. Comparative transcriptomics indicated no strong evidence of expression divergence in telomere maintenance genes for M. lewisii, suggesting rapid evolution shaped TERT protein sequence. In vivo activity of M. lewisii telomerase was investigated by analyzing F1 telomeres generated by crossing M. lewisii and M. verbenaceus, which doesnt have a functioning TR2. Results showed M. verbenaceus chromosome ends in the F1 had converted into M. lewisii telomeres, suggesting dominance of the M. lewisii telomerase. We demonstrate TERT-TR coevolution can have significant consequences on the evolution of plant telomeres. Significance statementTelomeres protect chromosome ends and are maintained by the telomerase complex. We discovered the catalytic component of the telomerase (TERT) was rapidly evolving in monkeyflowers (Mimulus) and studied the molecular consequences. In M. lewisii, TERT evolved lineage-specific amino acids to bind two sequence divergent telomerase RNA paralogs. Telomerase activity assay showed M. lewisii synthesized more telomere repeats compared to its sister species without the TR duplication, and transcriptomics indicated this was not due to a change in telomere maintenance gene expression. Genetic experiments in interspecies hybrids showed M. lewisii telomerase could convert chromosome ends in sister species into M. lewisii-like telomeres suggesting functional dominance. We show rapid evolution of the telomerase can have significant effects on telomere evolution.

Rosa gallica L., the French rose, is a perennial, tetraploid, heterozygous species that naturally propagates by seed and sucker. It occurs in the wild, primarily in Europe, and also exists as cultivated varieties. R. gallica cultivars were extensively bred and cultivated in France at the beginning of the 19th century. Although several hypotheses have been proposed regarding the species expansion based on historical records, none have been assessed using molecular data. Indeed, its genetic diversity has so far been investigated only at local or regional scales, hindering the identification of the evolutionary factors shaping its present-day distribution. Using 29 sequenced microsatellites, we genotyped a comprehensive sample of 1618 individuals, including wild R. gallica from 219 sites across the species range, rose cultivars, and specimens from other Rosa species. We then detected clonal lineages and characterized the range-wide genetic diversity and structure, aiming to disentangle the roles of natural and human factors in shaping the distribution of R. gallica, with particular focus on France. French diversity appears particularly structured compared to the rest of the range, suggesting multiple origins within France. Populations in South Alps, Central Eastern Europe, and Eastern France appear to have recolonized naturally from a single southern glacial refugium. In contrast, populations in the western part of France likely resulted from more recent natural or human-mediated dispersal. Finally, clonal lineages containing both wild and cultivated individuals were predominantly found in France, highlighting the role of human-mediated dispersal in 28 of the 98 French sites studied. These findings show that the present-day natural range of R. gallica was shaped primarily by post-glacial recolonization, but also reveal a contribution of human activities to its recent dispersal, particularly in France, where cultivated varieties were intensively bred and exchanged.

Frugivory is essential for the maintenance of tropical forests, influencing seed movement, spatial structure, and ecosystem functioning. Although frugivory interactions are well documented at local scales, we still know relatively little about how network structure varies across biogeographic regions, particularly in equatorial lowland forests. In this study, we compare frugivory networks in lowland forests in the Guiana Shield (South America), Malaysia (Southeast Asia), and Central Africa (Gabon and Cameroon). Because these sites occur at similar latitudes, they provide a useful setting to examine how evolutionary history and biogeographic context, rather than climate alone, might shape the structure of the network. We analyze patterns across regions, focusing on implications for frugivory interactions and, ultimately, seed dispersal. Guided by the Integrative Hypothesis of Specialization (IHS), we focus on the role of megafauna and other large-bodied frugivores, which are important consumers of large fruits and dispersers of large seeds, and have experienced different historical trajectories in each continent. We also considered whether smaller frugivores, such as rodents, might partially compensate for the loss of Megafauna in the Americas through seed handling and caching behavior. The pantropical comparison revealed structural and functional differences between the studied frugivory networks. In the Guianas, lacking megafauna, the frugivory network showed strong modularity with nestedness within modules, forming a marked compound topology. Alternatively, in Central Africa, with elephants and great apes, the network showed intermediate modularity and nestedness, reflecting a balance between local compartmentalization and regional integration. Finally, in Malaysia, the network was predominantly nested, with a few highly connected species, especially figs and flying foxes, that linked most of the partners. Across continents, frugivory networks tend toward a compound topology, but the relative influence of modular and nested components shifts with biogeography, evolution, and ecology.

This study aims to quantify the genetic similarity of different species (from fish to humans) to the human reference genome (pp6, Homo sapiens.GRCh38) based on the allele presence/absence patterns of 33 language/cognition related gene SNV loci, identify key breakpoints during evolution, and evaluate the enrichment of language and cognition genes at these breakpoints. We designed a similarity calculation method relying on binary features (four columns for A/T/C/G), adopted five difference/distance measures (Sorensen, Rogers, Nei, Reynolds, and Hellinger), and converted them into similarity values (1/(1+distance)). For each method, samples were independently ranked, the first derivative of similarity was computed, and the top 12 peaks were selected as candidate breakpoints. Results show that the similarity curves from the five methods are highly consistent (correlation coefficients >0.9), with major peaks concentrated at positions 355, 363, 381, 382, 390, 400, etc., where the corresponding samples are predominantly ancient hominins and primates. Furthermore, we defined 13 peak groups (starting positions 355-401). For each peak within a group, pairwise SNV differences between the peak apex sample and its immediate left neighbor were compared, and the intersection F_INTERSECTION (shared differential loci) was obtained. For each F_INTERSECTION, we calculated the proportions of language genes and cognition genes. In addition, we computed the differential sets between adjacent groups' F_INTERSECTION to trace the gradual emergence of new loci. In F_INTERSECTION, language genes accounted for an average of 59.5%, and cognition genes for an average of 62.9%. The proportion of language genes reached a peak at position 383 (61.2%), while cognition genes peaked at position 386 (64.9%). High frequency peak samples include c25, c27, and ja2, suggesting that language cognition genes may have undergone independent intensification during Eurasian evolution. Differential analysis between adjacent F_INTERSECTION revealed a stepwise acquisition of new loci from position 355 to 401, with three bursts of newly added loci along the entire evolutionary axis. This study provides a quantitative framework based on similarity curves, offers a novel molecular perspective for understanding the evolution of language and cognitive abilities, and highlights the potential importance of East Asian archaic hominins in the evolution of language cognition genes.

1. Plants produce a great chemodiversity, which is the diversity of specialized metabolites (SMs). These SMs are produced in complex metabolic pathways and play an important role in inter-species interactions. There are numerous hypotheses about the evolutionary processes which brought about and maintain chemodiversity. Some have been partially tested in lab and field studies. However, some of their assumptions and predictions are better tested by quantitative modeling, and so far no quantitative model has investigated the role of metabolic pathways. 2. To close this gap, we developed an individual-based model for metabolic pathway evolution. It models enzymes creating metabolites with various modifications. Enzymes undergo inheritance and mutation. We used the model to compare the screening and interaction diversity hypotheses. 3. The screening hypothesis predicts promiscuous enzymes, genetic drift, the presence of many non-beneficial metabolites, and high metabolite richness. The interaction diversity hypothesis predicts specialized enzymes, selection, the almost exclusive presence of beneficial metabolites, and situation- dependent metabolite richness. We found that the patterns predicted by the screening hypothesis did not occur, while those predicted by the interaction diversity hypothesis did. 4. This provides reason to favor the interaction diversity hypothesis over the screening hypothesis when connecting empirical results to their evolutionary context

Motivation: The Interactive Tree of Life (iTOL) is widely used to display and annotate phylogenetic trees, but managing its format-sensitive annotation files impede reproducible high-throughput analyses. Among the maintained Python packages and versions evaluated, none combined template generation, taxonomic monophyly assessment and iTOL batch operations. Results: PyiTOL validates inputs, generates 31 iTOL template schemas (22 accepted by the live batch uploader), performs LCA-based monophyly classification with nested-monophyly detection, sampling-completeness states and polyphyletic subgroup decomposition, plus API upload and session replay. On a topology-constructed benchmark, all calls matched prespecified labels for 4,389 groups; on a 700-genome tree, binary mono/non-mono calls agreed with ETE4 for 409 genera; 17,294 GTDB R232 genera were processed in about 17 s. Availability and Implementation: PyiTOL 1.0.3 (Python [≥]3.10; Linux, macOS and Windows) is MIT-licensed at https://github.com/ZengZichao/PyiTOL and archived with test data at Zenodo (https://doi.org/10.5281/zenodo.22106806).

Background: Reproducible taxonomic collapsing and geological-timescale annotation of time-calibrated phylogenetic trees in R often require coordination among several packages and repeated code for label parsing, clade validation, plotting, and export. Workflow-managed analyses additionally benefit from non-interactive configuration, predictable diagnostics, and machine-readable exit status. Results: We present Rclade, an R package that consolidates the multi-package coordination required for taxonomic collapsing into a streamlined, single-function interface. Rclade provides (1) custom ggproto objects (GeomPolygonStraight/GeomSegmentStraight) that bypass coord_munch() interpolation to achieve straight-edge rendering of collapsed triangles in circular layouts; (2) automatic detection and parsing of four taxonomic-label formats (GTDB, Silva, NCBI, embedded) plus user-supplied custom regex, with explicit input-validation contracts and parsing-accuracy evaluation on real and derived test sets; and (3) workflow embeddability through YAML configuration, library-mode APIs, and standard Unix exit codes. Benchmarks on synthetic and real datasets (200-10,000 synthetic tips and real reference trees up to 10,122 tips; 5 replicates at every scale under a unified fully rendered measurement protocol) show that the full-pipeline overhead is modest for interactive use (median {approx}0.87 s in-session rendering and {approx}8.4 s process-level wall-clock at 10,000 tips). Conclusions: Rclade is a convenience layer over the ggtree/deeptime ecosystem that reduces boilerplate while adding targeted technical improvements for circular-layout rendering and format heterogeneity management.

How morphological complexity evolves is still enigmatic. While there is evidence in algae and plants as well as animals that diversification of the repertoire of transcription factors (TF) is causative for evolution of organismal complexity, there are many examples from lineages that follow their own way of complexity evolution, for example by expansion of particular families. For land plants, correlation of the size of the TF complement with number of cell types (as a proxy for morphological complexity) has been shown, and several families were identified as candidates to drive complexity evolution. Here, we expand a previously available dataset of cell type numbers from 12 to 82 proteomes and introduce a four class body plan scheme. We find that the total TF complement correlates with the number of cell types of Archaeplastida (primary plastid bearing plants and algae). We used TabPFN (Tabular Prior-data Fitted Network) for binary (uni- vs. multicellularity) as well as for four class Bauplan classification. TabPFN is able to predict the morphological complexity with high accuracy. This approach allows to determine organismal complexity based on the gene space of an organism. Based on our results, we can confirm that plant morphological evolution is driven by gain and expansion of TF families.

Understanding the spatial patterns of biodiversity and their environmental determinants is fundamental to ecology and conservation, particularly in mountainous transitional zones where species assemblages can shift rapidly. Amphibians, as highly sensitive organisms, are excellent indicators of environmental change, yet quantitative assessments of their diversity in the Huangshan Mountains-a biodiversity hotspot at the junction of the Palaearctic and Oriental realms in southern Anhui, China-remain limited. We conducted systematic line-transect surveys across 200 grids (5 km x 5 km) during the spring and autumn of 2023 and 2024, recording a total of 10,342 individuals belonging to 23 species, 20 genera, 9 families, and 2 orders. Three species (Fejervarya multistriata, Microhyla fissipes, and Bufo gargarizans) were identified as dominant, and two nationally protected species (Hoplobatrachus chinensis and Andrias davidianus) were detected. Inter-annual alpha diversity did not differ significantly, but pronounced seasonal variation was observed, with spring supporting higher Shannon-Wiener and Pielou evenness than autumn in both years. Using generalized additive models and model averaging, we found that annual precipitation and the normalized difference vegetation index (NDVI) were consistently the strongest positive predictors of Shannon-Wiener diversity, Simpson dominance, and species richness, while species richness also declined significantly with increasing human footprint. Total beta diversity (Soerensen dissimilarity) was very high (0.981) and overwhelmingly driven by species turnover (97.86%) rather than nestedness. Partial Mantel tests and distance-based redundancy analysis further revealed that environmental distances-particularly annual precipitation-significantly shaped overall beta diversity and its turnover component after accounting for geographic distance. These results highlight the predominant role of climatic and vegetation gradients in structuring amphibian assemblages in this subtropical mountainous region, providing baseline data to inform local conservation strategies.

Mulberry trees have high economic and ecological value, and a robust molecular marker system plus germplasm genetic diversity analysis is critical for innovative utilization of high-quality medicinal and economic mulberry germplasm. Here, 51 mulberry samples were used to develop SNP primers via genome resequencing, with the SNP-PCR system optimized by single-factor and orthogonal assays. The phenotypic diversity and SNP molecular marker genetic diversity of 1-deoxynojirimycin (1-DNJ) in mulberry leaves were analyzed respectively, and the genetic correlation between molecular markers and phenotypic traits was evaluated by Mantel test. Tested germplasm showed marked 1-DNJ variation (0.4805-2.5300 mg/g, CV=0.4241), reflecting rich genetic diversity. The optimal SNP-PCR system included Buffer (containing Mg{superscript 2}+) 2.2 L, 2.5 mM dNTP 0.4 L, forward and reverse primers (10 mol{middle dot}L-1) totaling 2.75 L, Taq DNA polymerase (5 U{middle dot}L-1) 0.3 L, DNA (50 ng{middle dot}L-1) 1.1 L, and ddH2O 13.65 L. 23 highly polymorphic ones amplified 91 loci (81 polymorphic, 89.10% polymorphism rate). Genetic diversity analysis showed that the average genetic distance was 0.3010, and the average expected heterozygosity (H) and Shannon information index (I) reached 0.4667 and 0.3104 respectively, indicating that the genetic differentiation among the tested mulberry germplasms was significant and the population had a moderate to upper level of genetic diversity. UPGMA clustering divided 51 germplasms into 6 major groups at a genetic similarity coefficient of about 0.7, while phenotypic clustering based on 1-DNJ content divided them into 2 major categories and 4 subcategories, with high 1-DNJ germplasm clustered independently. Mantel correlation analysis showed that 6 SNP sites were significantly weakly correlated with 1-DNJ content (r < 0.3, p < 0.05), and can be used as candidate molecular markers for subsequent genetic analysis of 1-DNJ content.This study established a stable mulberry SNP-PCR system, Analyze the molecular genetic characteristics of mulberry germplasm and DNJ phenotypic variation rules respectively, and provide basic data for cluster comparison. and provided a scientific basis for marker database improvement, germplasm identification and molecular-assisted breeding.

Reconstructing ancestral morphologies on a phylogenetic tree is a central task in evolutionary morphometrics. Established reconstruction methods, including multivariate Brownian-motion approaches, rely on linear assumptions and do not directly model the correlations between landmarks within a shape, which can oversimplify the reconstructed morphology. The DICAROS method (Diffeomorphic Independent Contrasts for Ancestral Reconstruction of Shapes; Severinsen et al., 2026) instead fuses sibling shapes along branches with large-deformation diffeomorphic (LDDMM) landmark dynamics that model these correlations, so that ancestors remain on the shape manifold. DICAROS was shown to outperform ordinary least-squares, Brownian-motion, and penalized-likelihood reconstruction, particularly on non-symmetric trees. The dicaros package repackages that pipeline as a documented, pip-installable tool that runs on arbitrary landmark datasets from a single command. It handles 2D and 3D landmarks, Newick and NEXUS trees, a choice of Euclidean or Frechet species means, optional anchor-based alignment, and tips backed by a single specimen, and it returns the reconstructed shapes for all nodes together with the tree relabelled at its internal nodes. We demonstrate dicaros on two new datasets: a 2D leaf dataset (217 species) and a 3D guenon skull dataset (22 species).

Polytypic species with large ranges may harbor unrecognized diversity because taxonomy ranks populations differing subtly in plumage as subspecies. The Ruddy Foliage-gleaner (Clibanornis rubiginosus) exemplifies this problem. It ranges from Mexico to Brazil, with 15 subspecies, and forms a non-monophyletic complex with two congeners, yet its songs had not been compared. We measured ten spectral and temporal variables on 104 recordings covering 14 of 15 subspecies. Bayesian linear mixed models showed three song groups: eight subspecies west of the Andes share a single-note song, whereas Amazonian and Guianan populations add a short introductory note and sing longer, lower-pitched songs. Within this group, watkinsorum sings the lowest-pitched and longest song and is phylogenetically closer to C. cinnamomeigula than to its Amazonian neighbors. The third group is C. cinnamomeigula alone, a white-eyed taxon in an otherwise dark-eyed group. Its high-pitched, vibrato song resembles none other in the genus. One-note and two-note songs differ in kind without intermediates, and every two-note taxon sequenced to date falls in one clade, so C. rubiginosus is paraphyletic. We recognize four species, C. rubiginosus sensu stricto, C. cinnamomeigula, C. watkinsorum, and C. obscurus. This raises Clibanornis from five species to eight and divides its only polytypic species.

Complex biomechanical innovations are often treated as discrete evolutionary breakthroughs, yet their diversification within large radiations remains poorly understood. Beetle leg joints provide a striking example: some weevils possess screw-like coxa-trochanteral articulations in which rotation and axial displacement are mechanically coupled, resembling engineered screw-and-nut mechanisms. Whether these joints represent isolated mechanical extremes, discrete adaptive types or part of a broader continuum of phenotypic variation has remained unknown. Here we combine synchrotron X-ray microtomography, landmark-free atlas-based morphometrics, quantitative functional morphology and phylogenetic comparative analyses to examine the mesocoxa-trochanteral joint in 68 specimens representing seven sampled family-level groups across early-diverging and derived weevil lineages. We show that screw joint evolution combines continuous variation in trochanteral shape with a restricted set of mechanically plausible joint-character combinations, rather than forming sharply separated morphological classes. True screw-and-nut joints are not confined to a distinct region of morphospace, indicating that overall form and mechanical configuration are not necessarily coupled. The occurrence of this configuration in the early diverging Caridae shows that it is not restricted to more derived families. Three-dimensional helix fitting revealed a mosaic geometry, with winding angle showing the clearest relationship with overall shape and joint architecture, whereas axial pitch varied largely independently of shape, size and lineage. Together, these patterns show that screw joint components diversified with different degrees of evolutionary integration. These results recast the weevil screw joint from a singular biomechanical curiosity into a diversified evolutionary system. They suggest that complex functional structures can evolve through the gradual recombination and differential persistence of structurally constrained and evolutionary flexible components, rather than through a single shift from simple to fully specialized designs.

Among Madagascar primates, the sportive lemurs (family Lepilemuridae) have seen their species diversity increase from eight in 2005 to 26 in 2009 mostly by applying the phylogenetic species concept to DNA barcode data. Despite the genus being speciose, only one case of sympatry is known from northern Madagascar, where two sportive lemur species described based on low mtDNA divergence, Lepilemur ankaranensis and Lepilemur milanoii, were found to co-occur at the center of their joint distribution range. Here, to clarify the taxonomy of these two species and examine their sympatry, we apply an integrative taxonomic framework to genomic and morphological data from 84 individuals of L. ankaranensis and L. milanoii, encompassing their entire distribution range and the forest of Analafiana, beyond their southernmost limit. Using clustering, multivariate, and isolation by distance analyses, we find no evidence of a sympatric zone and show that despite clear genetic differentiation between regions, the genomic and morphological diversity of the L. ankaranensis, L. milanoii-Analafiana group is clinal and explained by geographic distance. These results clarify that L. milanoii is a junior synonym of L. ankaranensis and that the Analafiana forest population belongs to L. ankaranensis, extending its distribution. It further implies that the 'sympatric' zone, the Andrafiamena forest, hosts conspecific individuals with slightly differentiated mtDNA backgrounds, rather than slightly differentiated sympatric species. Lastly, we re-evaluate the IUCN conservation metrics of L. ankaranensis, which continue to qualify as Endangered (EN) under the B1ab(i-v) criteria.

Background and AimsDistylous polymorphisms (distyly) are adaptations for plants to improve efficiency of cross-pollination and reduce pollen wastage. Theoretical models suggest that distyly evolves from stylar monomorphism via an intermediate stage of stigma-height dimorphism (SHD), however, this evolutionary scenario is only observed in few distylous lineages. The genus Jasminum have many species exhibit either distyly or SHD, providing an ideal opportunity to test models of the evolution of distyly. Focusing on the evolution of distyly, we investigated floral morphs and evaluated the occurrence of distyly and SHD in Jasminum using phylogenetic reconstruction and morphological analysis. MethodsWe investigated floral morphs and evaluated the occurrence of distyly and stigma-height dimorphism through morphology observation. To perform phylogenetic analysis, we sequenced plastomes of forty species from Jasminum and Chrysojasminum using Illumina next generation sequencing, and then constructed phylogenetic tree using maximum likelihood method and Bayesian inference. Based on the phylogenetic tree, we inferred the ancestral floral type through ancestral reconstruction. Key ResultsOur results suggest that the distyly originated in the common ancestor of Jasminum and Chrysojasminum. SHD occurs in several Jasminum species scattered in different clades of the phylogenetic tree, suggesting multiple independent reversions from distyly back to SHD. However, this transition is not associated with the loss of ancillary polymorphisms, as all the species examined in this study well retain dimorphic traits for other floral organs such as pollens and stigmas. Moreover, most species maintain strict heteromorphic self-incompatibility, while J. officinale has lost or at least partially lost self-incompatibility, suggesting that distyly is not always linked to self-incompatibility. ConclusionsIn Jasminum, breakdown of distyly resulted in evolutionary transitions to stigma-height dimorphism for multiple times, suggesting that distyly is not a stable floral polymorphism under certain selective forces. These findings advance our understanding on the evolution of distyly and plant reproductive systems.

Leaf shape is a fundamental trait of plant ecological strategies, influencing biotic interactions and ecosystem functioning. However, established quantitative metrics fail to capture subtle variations and irregularities, require user-based reference points or are challenging to compare among taxa with broadly different leaf shapes. In addition, established metrics typically conflate (aggregate) leaf edge complexity and macro-shape complexity, despite their independent functional significance and genetic foundations. Here, we introduce an entropy-based framework to quantify two new complexity metrics: edge complexity and macro-shape complexity. Based on three case studies, we show that these metrics outperform aggregate metrics in predicting Quercus robur chemical traits, provide more intuitive interspecific classifications, and strongly align with human perception. In addition, edge and macro-shape complexity show high complementarity, while aggregate metrics are highly redundant and typically strongly related to leaf area. Emerging as the strongest predictor of leaf chemistry and key visual cue for complexity as perceived by humans, the effects of edge complexity highlight the under-appreciated functional significance of leaf margins. Our framework and the proposed entropy-based complexity metrics thus promise to help unlock the potential of growing digital image archives of leaves, including images from herbaria and fossils, and are technically readily applicable to shapes of algae, bacteria, pollen, and beyond. The accompanying package ShapeComplexity enables the broad application of entropy-based metrics, providing a powerful tool to explore how the shape of organisms and biological structures influences ecological strategies, biotic interactions, and ecosystem functioning while tracking spatial and temporal variation.